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Abstract. A hot central star illuminating the surrounding ionized H 
II region usually produces very rich atomic spectra resulting from basic 



(N 



i-j^ ' atomic processes: photoionization, electron-ion recombination, bound- 

bound radiative transitions, and collisional excitation of ions. Precise 



OS 



diagnostics of nebular spectra depend on accurate atomic parameters 
for these processes. Latest developments in theoretical computations 
are described, especially under two international collaborations known 
as the Opacity Project (OP) and the Iron Project (IP), that have yielded 
' accurate and large-scale data for photoionization cross sections, transi- 

tion probabilities, and collision strengths for electron impact excitation of 
most astrophysically abundant ions. As an extension of the two projects, 
£vq ■ a self-consistent and unified theoretical treatment of photoionization and 

CO [ electron-ion recombination has been developed where both the radiative 

. and the dielectronic recombination processes are considered in an uni- 

fied manner. Results from the Ohio State atomic-astrophysics group, 
^ ■ and from the OP and IP collaborations, are presented. A description 

of the electronic web-interactive database, TIPTOPBASE, with the OP 
O ' and the IP data, and a compilation of recommended data for effective 

collision strengths, is given. 
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c3 



1. Introduction 

A planetary nebula may be thought of as an 'astrophysical laboratory' of atomic 
emission-line spectra. The spectra enable determination of its temperature, 
density, and abundances of elements. The precise spectral analyis require accu- 
rate atomic parameters for the radiative are collisional processes in the nebular 
plasma. The four basic atomic processes that dominant the nebular plasma are: 

i) Radiative bound-bound transition for excitation or de-excitation: 

X +z + hv # X +z * , 

where the ion, X, is of charge, Z. 

ii) Photoionization (PI) by absorption of a photon: 

X +z + hv^ X +z+1 + e, 

The inverse process is the electron- ion radiative recombination (RR). 

iii) Autoionization (AI) and dielectronic recombination (DR): 
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e + X +z Al 
X +z - l + hv DR 



where the intermediate doubly excited state, autoionizing state, introduces res- 
onances in the atomic processes. The inverse process of DR is photoionization 
via the autoionizing state. 

iv) The collisional process of electron-impact excitation (EIE) 



is one of the primary processes for spectral formation in astrophysical plasmas. 

Each process needs to be treated separately to obtain the relevant atomic 
parameters. Sample results for each process, obtained mainly under the OP 
(1995, 1996) and the IP (1993) for accurate atomic data and stellar opacities 
will be presented later. Astrophysical model applications, such as for plasma 
opacities, need a huge amount of atomic data as these deal with a large number 
of atomic levels over wide energy ranges. Consistent sets of atomic parameters 
can be obtained if the atom or the ion is described by the same wavefunction 
expansion in each process. This also reduces the uncertainty in applications 
involving different processes and approximations. The Close Coupling (CC) R- 
matrix mehtodology employed by the OP and IP enable the computation of such 
self-consistent sets of atomic parameters. 

The large amount of atomic data and opacities obtained under the OP 
and the IP are available electronically through the existing database, TOPbase 
(Cunto et al. 1993), and through its planned extension TIPTOPBASE (C. 
Mendoza and the OP/IP team). 

2. Theory 

All calculations for the various atomic parameters of the dominant atomic pro- 
cesses are carried out using the accurate and powerful R-matrix method in the 
close-coupling approximation (e.g. Burke & Robb 1975, Seaton 1987, Berring- 
ton et al. 1987, Berrington et al. 1995). The total wavefunction for a (N+l) 
electron system in the CC approximation is described as: 



where Xi is the target ion or core wavefunction in a specific state SiL^i or level 
Jj7Tj, 9i is the wavefunction of the interacting (N+l)th electron in a channel 
labeled as SiLi(Ji)iTi kf£i(SLir or Jir); kf is the incident kinetic energy. is 
the correlation functions of (e+ion) system that compensates the orthogonality 
condition and short range correlation interations. The complex resonant struc- 
tures in photoionization, recombination, and in electron impact excitation are 
included through channel couplings. 

Relativistic effects are included through Breit-Pauli approximation in inter- 
mediate coupling. The (N+l)-electron Hamiltonian in the Breit-Pauli approxi- 



e + X +z - e' + 



X +z ** 



N 

^e(g + ion) = A*^2 Xi{i- on )Qi + c 3^j( e + ^ on ) 
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mation, as adopted in the Iron Project, is 



H 



BP 

N+l 



ttNR i u mass i zjDar i zjso 
H N+l + -"/V-4-1 + -"/V+1 + tl 



l N+l 



where non-relativisitc Hamiltonian is 
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NR 
N+l 
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JV+li 
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is the mass correction term, iS^+i * s the Darwin term, Hj^ +1 is the 
spin-orbit interaction term. Spin-orbit interaction splits the LS terms into fine- 
structure levels labeled by Jtt where J is the total angular momentum. Solutions 
of the Schrodinger equation, H^+i& = which becomes a set of coupled equa- 
tions with the CC expansion, give the bound wavefunctions, *$>b, for negative 
energies, E < 0, and continuun wavefunction, typ, for positive energies, E > 0. 
The transition matrix elements for various atomic processes are: 

< ^^IIDII^f > for photoionization and recombination, 

< \I/£||D||\I/b' > for oscillator strength, 

< typ\H(e + ion)\^pi > for electron impact excitation, 
where D is the dipole operator. 

The transition matrix element with the dipole operator can be reduced to 
the generalized line strength defined, in either length or velocity form, as 



I N+l 
\ 3=1 



Sy = CO 



-2 



JV+1 







(4) 



where u is the incident photon energy in Rydberg units, and <fj and \Pf are the 
wave functions representing the initial and final states, respectively. 

The oscillator strength and the transition probability Aji for the bound- 
bound transition can be obtained from S as 



fij — g^ - ^' Aji{a.U. 



2° g 3 ^ 



(5) 



where Eji is the transition energy, a is the fine structure constant, and gi, §j 
are the statistical weight factors of the initial and final states. The lifetime, 
Tj, of a level j can be obtained from the transition probabilities decaying to 



the lower levels as Tj = (J2i^-ji{ s 1 )) X > where Aji(s r ) 



Ajija.u.) 



TO 



and tq = 



2.4191 x l(T 17 s. 

The photoionization cross section (a pi) is proportional to the generalized 
line strength (5), 

47Fl Q 

api = loS. (6) 

3c gi 

The calculations of total electron-ion recombination employs a new unified 
treatment (Nahar & Pradhan 1994,1995). The method considers the infinite 
number of recombining states and incorporates the radiative and dielectronic 
recombinations, RR and DR, in a unified manner. The contributions to re- 
combination from states with n < 10 are obtained from opi using principle of 
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detailed balance, while the resonant contributions from states with 10 < n < oo 
are obtained from an extension of the DR theory of Bell & Seaton (1985). 

The recombination cross section, (Jrc, is related to opi through the princi- 
ple of detailed balance, 

9i _ffuP_ 

<*RC = OPI—-TT, n 9 n - {<) 

gj 4:TT z m z c z v z 

The recombination rate coefficient, oirc, is obtained as 

aRc(T) = / vf(v)a RC dv, (8) 
Jo 

where f(v) is the Maxwellian velocity distribution function. The total o.rc is 
obtained from contributions from infinite number of recombined states. 

The collision strength for transition by electron impact excitation from the 
initial state of the target ion SiLi to the final state SjLj is given by 

n(SiLi - SjLj) = ^J2 E( 25 + !)( 2L + l)\S SL *(SiLiU ~ SjLjlj)?, (9) 

SLn lilj 

where S is the scattering matrix. The effective collision strength or the Maxwellian 
averaged collision strength can be obtained as 

T(T)= rn^e^diej/kT). (10) 
Jo 

which give the excitation rate coefficient, qij{T) = (8.63x /ujiT l l 2 )e- E ^ /kT T{T) 
in cm 3 s _1 . where T is in K, Eij = Ej — Ei, E^ < Ej are in Rydbergs (1/kT = 
157885/T), and j is the excited upper state. 



3. Results and Discussions 

Sample results for each atomic process are presented, e.g., for /-values, api, 
otR(T), and T(T), in the following subsections. 

3.1. Bound-Bound transitions 

Due to fine structure, extensive sets of data for /- and ^4-values can be obtained 
for a large number of bound-bound transitions using the Breit-Pauli R-matrix 
method. In contrast to dipole allowed LS multiplets, BPRM method includes 
relativistic effects, and consideration of both the dipole allowed and intercombi- 
nation transitions. The fine structure energy levels are analysed with quantum 
defects to obtain the spectroscopic identification, (Ct St L t Jt irtn£)SLJ ir where 
(Ct St L t J t Tit) denotes the core configuration, spin, orbital, total angular mo- 
menta and parity, nl is configuration of the outer or valence electron, SL are 
the possible spin and orbital angular momenta, and Jit are the total angular 
momentum and parity of the (N+l)-electron system. The fobidden quadrupole 
(E2) and magnetic dipole (Ml) transitions are treated with atomic structure 
calculations using codes such as SUPERSTUCTURE (Eissner et al. 1974). 
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Accurate fine structure transition probabilities have been obtained for a 
number of ions such as: 

Li-like: C IV, N V, VI, F VII, Ne VIII, Na IX, Mg X, Al XI, Si XII, S XIV, 

Ar XVI, Ca XVIII, Ti XX, Cr XXII, Ni XXVI 

Fe ions: Fe V, Fe XVII, Fe XXI, Fe XXIII, Fe XXIV, Fe XXV 

Other ions: C II, C III, IV, S II, Ar XIII 

The first large-scale application of the BPRM method for a complex ion, 
Fe V, resulted in 3865 fine structure energy levels, compared to 179 observed, and 
about 1.46 million /- and A— values for dipole allowed and intercombination (El) 
transitions and 362 forbidden (E2,M1) transitions (Nahar et al. 2000). Table 
I presents a sample with complete spectroscopic identification of the transition 
array 3d 4 ( 4 D) - 3d 3 ( 4 F)4p{ 5 F°) in Fe V. 



Table I: Transition probabilities of Fe V. g=2J+l. 



leveL 



lev el j 



g { gj Ei(Ry) Ej(Ry) 



f 



A(sec~ l ) 



3d 4 ( 5 D) 
3d 4 ( 5 D) 
3d 4 ( 5 D) 
3d 4 ( 5 D) 
3d 4 ( 5 D) 
3d 4 ( 5 D) 
3d 4 ( 5 D) 
3d 4 ( 5 D) 
3d 4 ( 5 D) 
3d 4 ( 5 D) 
3d 4 ( 5 D) 
3d 4 ( 5 D) 



3d 3 ( 4 F)Ap( 5 F°) 1 3 5.5132 3.1644 2.154E-01 3.18E+09 

3d 3 ( 4 F)4p( 5 F°) 3 3 5.5119 3.1644 3.790E-04 1.68E+07 

3d 3 ( 4 F)4p( 5 F°) 3 5 5.5119 3.1496 1.358E-03 3.65E+07 

3d 3 ( 4 F)4p( 5 F°) 5 3 5.5094 3.1644 4.617E-02 3.40E+09 

3d 3 ( 4 F)4p( 5 F°) 5 5 5.5094 3.1496 5.967E-02 2.67E+09 

3d 3 ( 4 F)4p( 5 F°) 5 7 5.5094 3.1443 1.462E-02 4.69E+08 

3d 3 ( 4 F)4p( 5 F°) 7 5 5.5058 3.1496 6.895E-03 4.30E+08 

3d 3 ( 4 F)4p( 5 F°) 7 7 5.5058 3.1443 5.889E-02 2.64E+09 

3d 3 ( 4 F)4p( 5 F°) 9 7 5.5015 3.1443 1.966E-03 1.13E+08 

3d 3 ( 4 F)Ap( 5 F°) 7 9 5.5058 3.1391 3.262E-02 1.14E+09 

3d 3 ( 4 F)Ap( 5 F°) 9 9 5.5015 3.1391 5.139E-02 2.30E+09 

3d 3 ( 4 F)Ap( b F°) 9 11 5.5015 3.1343 7.548E-02 2.78E+09 



3d 4 ( 5 D) 3d 3 { 4 F)4,p{ 5 F°) 25 35 5.5055 3.1451 1.068E-01 3.42E+09 



3.2. Photoionization and Recombination 

The photoionization cross sections, dpi, are calculated including autoionizing 
resonances that can enhance the background cross sections considerably. Fig. 1 
shows the photoionization cross sections of the gound states of Fe I to Fe V 
(Bautista & Pradhan 1998). Extensive resonances dominate the cross sections 
for these complex ions. The enhancement in the backgound is up to three orders 
of magnitude for Fe I, over an order or magnitude for Fe II, and ~ 50% for Fe 
III. 

For simpler systems the accuracy of theoretical cross sections can be tested 
against experimental data. Very precise measurements of a pi are now being 
carried out at a few places, such as the Advanced Light Source at Berkeley (e.g. 
Covington et al. 2001), the merged photon-ion beam set-up at Aahus University 
(Kjeldsen et al. 1999), and a synchroton based experiment at University of Paris- 
Sud ( Wuilleumier et al., private communication). An example of the comparison 
of ground state photoionization cross sections of C II, with very good agreement 
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Figure 1. Photoionization cross sections, crp/, of the ground states 
of Fe I - Fe V. 




(L>) Theory: BPRM 
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Figure 2. Comparison of theoretical and experimental photoion- 
ization cross sections, crpi, of the ground state fine structure levels 
2s 2 2p( 2 P 1 ° /23/2 ) of C II (Nahar 2002). 



with experiment, is shown in Fig. 2. Theoretical calculations include relativistic 
fine structure. 

An example of the total unified recombination rate coefficient a/?(T) for 
O III is presented in Fig. 3 (solid) (Nahar 1998). The earlier results are RR rates 
(dashed, Pequignot et al. 1991), low temperature DR rates (dotted, Nussbaumer 
& Storey 1983), high temperature DR rates by Badnell & Pindzolla (1990, short- 
dash-long-dash) and by Shull & Steenberg (1982, dot-desh). Differences can be 
noticed between the present unified values and the sum of the earlier (RR+DR) 
results, especially in the low and intermediate temperature. Compared to the 
earlier Shull & Steenberg results, autoionization into excited levels in the high 
temperature region reduces the recombination rates significantly. 

The unified treatment for the total recombination provides rates that are 
valid over a wide range of temperatures for all practical purposes in contrast to 
addition of RR and DR rates obtained using different approximations for dif- 
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Figure 3. Comparison of the unified total recombination rate coeffi- 
cients (ojr(T) (solid) with the previous calculations. 



ferent temperature ranges. The method enables obtaining self-consistent sets of 
photoionization and recombination cross sections by using the identical wave- 
function expansion for both the inverse processes. Following is the list of atoms 

and ions for which self-consistent sets of a pi and an{T) have so far been obtained 

for over 45 ions (e.g. Nahar & Pradhan 1997, http : // www.astronomy.ohio — state.edu/ ^p radhan): 

Carbon: C I, C II, C III, C IV, C V, C VI 

Nitrogen: N I, N II, N II, N IV, N V, N VI, N VI 

Oxygen: O I, O II, O III, O IV, O V, O VI, O VII, O VII 

Iron: Fe I, Fe II, Fe III, Fe IV, Fe V, Fe XIII, Fe XVII, Fe XXI, Fe XXIV, Fe 

XXV, Fe XXVI 

C-like: F IV, Ne V, Na VI, Mg VII, Al VIII, Si IX, S XI, Ar XIII, Ca XV 
Other ions: Si I, Si II, S II, S III, Ar V, Ca VII, Ni II 

The data include the state specific recombination rates for hundreds of bound 
levels with n< 10 for each ion. 

The self-consistent sets of photoionization/recombination data include new 
photoionization cross sections that are generally an improvement over the OP 
data since more extensive and accurate eigenfunction expansions are employed. 



3.3. Electron Impact Excitation 

Reviews and compilations of available theoretical data sources are: An evaluated 
compilation of theoretical data sources for electron-impact excitation of atomic 
ions, (A.K. Pradhan & J.W. Gallagher, Atomic Data and Nuclear Data Tables, 
52, 227, 1992), and Electron Collisions with Atomic Ions (A.K. Pradhan & H.L. 
Zhang, in LANDOLT-BORNSTEIN Volume Atomic Collisions, Ed. Y. Itikawa, 
Springer- Verlag, in press). 

A table of recommended data for effective collision strengths and A- values 
for nebular ions is available on-line from www.astronomy.ohio-state.edu/~pradhan 

The aim of the IP is to compute collisional data for the iron-peak elements 
in various ionization stages. Of the 50 publications in the "Atomic data from 
the IRON Project" series in Astronomy and Astrophysics journal, most are on 
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Figure 4. Collision strength for the first transition in Fe VI 
r2( 4 i ? 3/ 2 — 4 -?5/2) from the Iron Project CC R-matrix calculation (Chen 
& Pradhan 1999); previous results (• - Nussbaumer & Storey 1978, x 
- Garstang et al. 1978) neglect resonances. 

collisional strengths. At present the Iron Project website is maintained by Keith 
Butler at: www.usm.uni — muenchen.de/people/ip/iron —project.html 

An example of IP calculations for iron ions is the recent work on Fe VI (Chen 
& Pradhan 1999). Fig. 4 presents the detailed collision strength VL{E) showing 
extensive resonance structure. Previous calculations neglecting resonances seri- 
ously underestimate the Maxwellian averaged effective collision strength T(T) 
(Eq. 10) by several factors. 



4. TIPTOPBASE: atomic radiative and collisional data 



The existing electronic ar chive for the OP data, TOPbase (|/i£tp : / /heasarc.gsfc.nasa.goi 



at Goddard, NASA and fittp : / /vizier. u — strasbg.fr/OP.html at CDS) con- 



tain (i) photoionization cross sections of bound LS terms for ~ 200 ions with Z 
= 1 - 14, 16, 18, 20, 26, (ii) energy levels and transition probabilities (~ 10 7 f- 
values), (iii) monochromatic and Rosseland mean opacities (at CDS only) 

The new database, TIPTOPbase (under development, C. Mendoza & the 
OP/IP team) will have radiative and collisional data from the OP and the IP. 
The data include: 

(i) All data from TOPbase (and replacement of improved data), (ii) col- 
lisional data for iron and iron peak elements, (iii) new elements through all 
ionization stages: P, CI, K, as well as selected ions such as Ni II and Ni III, (iv) 
updated sets of radiative data from new CC calculations, (v) "Tail" photoion- 
ization cross sections at high energies that include inner-shell ionization, (vi) 
total and level specific recombination rate coefficients, (vi) additional data for 
/-values including inner-shell excitations in iron ions Fe VIII - XIII ("PLUS" 
data) calculated with SUPERSTRUCTURE, (vii) radiative data (crpj and /- 
values) for fine structure levels including relativistic effects, (viii) on-line com- 
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putational facilities for opacities and radiative accelerations for user-specified 
mixtures of elements ("customized" opacities and radiative forces) (see IP in 
Extended Abstracts). 

5. Conclusion 

The current status of large-scale ab initio close coupling R-matrix calculations for 
radiative and collisional processes is reported. The Iron Project Breit Pauli R- 
matrix radiative calculations include large numbers of dipole allowed and inter- 
combination transitions. Self-consistent sets of atomic data for photoionization 
and unified (electron-ion) recombination (including RR and DR) are obtained, 
and should yield more accurate photoionization models. Work is in progress for 
heavy ions of the iron group elements. 

Acknowledgments. This work is supported partially by the U.S. National 
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